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We present recent works [15,16] on the thermodynamic limit of quantum 
Coulomb systems. We provide a general method which allows to show the 
existence of the limit for many different systems. 

Ordinary matter is composed of electrons (negatively charged) and nu- 
clei (positively charged) interacting via Coulomb forces. The potential be- 
tween two particles of charges z and z' located at x and x' in R 3 is 



\x — x \ 

There are two difficulties which occur when trying to describe systems com- 
posed of electrons and nuclei. Both have to do with the physical problem 
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of stability of quantum systems. 

The first is due to the singularity of l/\x\ at 0: it is necessary to explain 
why a particle will not rush to a particle of the opposite charge. One of the 
first major triumphs of the theory of quantum mechanics is the explanation 
it gives of the stability of the hydrogen atom (and the complete description 
of its spectrum) and of other microscopic quantum Coulomb systems, via 
the uncertainty principle. Stability means that the total energy of the con- 
sidered system cannot be arbitrarily negative. If there was no such lower 
bound to the energy it would be possible in principle to extract an infinite 
amount of energy. One often refers to this kind of stability as stability of 
the first kind [19,20]. If we denote by E(N) the ground state energy of the 
system under consideration, for N particles stability of the first kind can 
be written 

E(N) > -co. (1) 

In proving (1) for Coulomb systems, a major role is played by the un- 
certainty principle which for nonrelativistic systems is mathematically ex- 
pressed by the critical Sobolev embedding H 1 (M. 3 ) L 6 (R 3 ). The latter 
allows to prove Kato's inequality 

V£>0, J_< e (_A) + -, 

\X\ £ 

which means that the Coulomb potential is controlled by the kinetic energy. 

The second issue concerns the slow decay of l/\x\ at infinity and this 
has to do with the macroscopic behavior of quantum Coulomb systems. It 
is indeed necessary to explain how a very large number of electrons and 
nuclei can stay bounded together to form macroscopic systems, although 
each particle interacts with a lot of other charged particles due to the 
long tail of the Coulomb interaction potential. Whereas the stability of 
atoms was an early triumph of quantum mechanics it, surprisingly, took 
nearly forty years before the question of stability of everyday macroscopic 
objects was even raised (see Fisher and Ruelle [11]). A rigorous answer 
to the question came shortly thereafter in what came to be known as the 
Theorem on Stability of Matter proved first by Dyson and Lenard [8] . 

The main question is how the lowest possible energy E(N) appearing 
in (1) depends on the (macroscopic) number N of particles in the object. 
More precisely, one is interested in proving a behavior of the form 

E(N) eN. (2) 

This behavior as the number of particles grows is mandatory to explain why 
matter does not collapse or explode in the thermodynamic limit. Assume 
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that (2) does not hold and that for instance E(N) <~;v->oo cN p with p ^ 
1. Then \E(2N) — 2E(N)\ becomes very large as N » 1. Depending on 
p and the sign of the constant c, a very large amount of energy will be 
cither released when two identical systems are put together, or necessary 
to assemble them. The constant e in (2) is the energy per particle. 

Stability of Matter is itself a necessary first step towards a proof of (2) 
as it can be expressed by the lower bound 

E(N) > -kN. (3) 

Put differently, the lowest possible energy calculated per particle cannot be 
arbitrarily negative as the number of particles increases. This is also often 
referred to as stability of the second kind [19,20]. 

A maybe more intuitive notion of stability would be to ask for the 
volume occupied by a macroscopic object (in its ground state). Usually 
this volume is proportional to the number of particles N. Denoting by O a 
domain in M 3 which is occupied by the system under consideration and by 
E(Q) its (lowest possible) energy, (2) then reads 

E(n) ~| nHoo e|fi| (4) 

where is the volume of ft. Stability of the second kind is expressed as 

E(n) > -k\SI\. (5) 

Instead of the ground state energy, one can similarly consider the free 
energy F(fl,/3,fj,) at temperature T —1/(3 and chemical potential /i. One 
is then interested in proving the equivalent of (4) 

F(fi,/?,p)~|nHoo/(MM (6) 

where /(/3,/i) is the free energy per unit volume. 

Large quantum Coulomb systems have been the object of an important 
investigation in the last decades and many techniques have been developed. 
A result like (3) (or equivalently (5)) was first proved for quantum electrons 
and nuclei by Dyson and Lcnard [8] . After the original proof by Dyson and 
Lcnard several other proofs were given. Licb and Thirring [27] in particular 
presented an elegant and simple proof relying on an uncertainty principle 
for fermions. The different techniques and results concerning stability of 
matter were reviewed in several articles [19-21,29,33]. 

It is very important that the negatively charged particles (the electrons) 
are fermions. It was discovered by Dyson [7] that the Pauli exclusion prin- 
ciple is essential for Coulomb systems: charged bosons are alone not stable 
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because their ground state energy satisfies E(N) ~ — C7V 7 / 5 , as was proved 
later [4,26,32]. 

A result like (2) (or equivalently (4)) was first proved by Lieb and 
Lebowitz [22] for a system containing electrons and nuclei both consid- 
ered as quantum particles, hence invariant by rotation. Later Fefferman 
gave a different proof [9] for the case where the nuclei are classical particles 
placed on a lattice, a system which is not invariant by rotation. 

In a recent work [15,16], we provide a new insight in the study of the 
thermodynamic limit of quantum systems, by giving a general proof of (4) 
or (6) which can be applied to many different quantum systems includ- 
ing those studied by Lieb and Lebowitz [22] or Fefferman [9], and others 
which were not considered before. Our goal was to identify the main gen- 
eral physical properties of the free energy which are sufficient to prove the 
existence of the thermodynamic limit. However, for the sake of simplicity 
we will essentially address the crystal case in this paper and we refer to our 
works [15,16] for a detailed study of the other cases. 

In proving the existence of the thermodynamic limit of Coulomb quan- 
tum systems, the most difficult task is to quantify screening. Screening 
means that matter is arranged in such a way that it is essentially locally 
neutral, hence the electrostatic potential created by any subsystem decays 
much faster than expected. This effect is the main reason of the stability 
of ordinary matter but it is very subtle in the framework of quantum me- 
chanics because the particles are by essence delocalized. In our approach, 
we shall heavily rely on an electrostatic inequality which was proved by 
Graf and Schenker [12,13] and which serves as a way to quantify screening. 
It was itself inspired by previous works of Conlon, Lieb and Yau [4,5], for 
systems interacting with the Yukawa potential. Fefferman used a similar 
idea in his study of the crystal case [9]. 

Like in previous works, our method consists in first showing the existence 
of the limit (6) for a specific domain A which is dilated (and possibly 
rotated and translated). Usually A is chosen to be a ball, a cube or a 
tetrahedron. In the applications [16] we always choose a tetrahedron as 
we shall use the Graf-Schcnkcr inequality [13] which holds for this type of 
domains. The second step consists in showing the existence of the limit (6) 
for any (reasonable) sequence of domains {£l n } such that |0„| — ► oo. This 
is important as in principle the limit could depend on the chosen sequence, 
a fact that we want to exclude for our systems. We shall specify later what 
a "reasonable" sequence is. Essentially some properties will be needed to 
ensure that boundary effects always stay negligible. 
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It is to be noticed that our method (relying on the Graf-Schcnkcr in- 
equality) is primarily devoted to the study of quantum systems interacting 
through Coulomb forces. It might be applicable to other interactions but 
we shall not address this question here. 

Proving a result like (4) or (6) is only a first step in the study of the 
thermodynamic limit of Coulomb quantum systems. An interesting open 
problem is to prove the convergence of states (or for instance of all fc-body 
density matrices) and not only of energy levels. For the crystal case, conver- 
gence of the charge density or of the first order density matrix was proved 
for simplified models from Density Functional Theory or from Hartree-Fock 
theory [3,25]. A result of this type was also proved for the Hartree-Fock ap- 
proximation of no-photon Quantum Electrodynamics [14]. 

Another (related) open question is to determine the next order in the 
asymptotics of the energy in the presence of local perturbations. Assume 
for instance that the crystal possesses a local defect modelled by a local 
potential V and denote the ground state energy in the domain by E v (Q) . 
Since V is local, it does not contribute to the energy in the first order of 
the thermodynamic limit. One is then interested in proving a behavior like 
E V (Q) = E°(Q) + f(V) + o(l)|n|- > oo- Such a result was recently proved 
for the reduced Hartree-Fock model of the crystal with the exchange term 
neglected [2]. This includes an identification of the function f(V). This 
program was also tackled for the Hartree-Fock model (with exchange term) 
of no-photon Quantum Electrodynamics [14]. 

The present paper is organized as follows. In the first section we intro- 
duce the model for the crystal and state our main theorem. In Section 2, 
we briefly describe two other quantum systems which we can treat using 
our method. Section 3 is devoted to the presentation of our new approach, 
in a quite general setting, together with hints on how it can be applied to 
the crystal case. 



1. The Crystal Case 

For simplicity, we put identical nuclei of charge +1 on each site of 1? . The 
results below can be generalized to any periodic system. Let be a bounded 
open set of M 3 and define the ./V-body Hamiltonian in il by 




i=i 
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where 

N 

w = E E j^+l E ^7^+^ E jjnhr 

Here — A is the Dirichlet Laplacian on (we could as well consider another 
boundary condition) . The Hamiltonian i?^ acts on iV-body fcrmionic wave- 
functions ^f(xi, .., xn) e /\^ L 2 (Q). Stability of the first kind states that 
the spectrum of Hg is bounded from below: 

E ° = = inf ff Af > -oo. 

ll*lll,2=l 

We may define the ground state energy in Q, by 

E(Sl) := inf Eg. (7) 

N>0 

It is more convenient to express (7) in a grand canonical formalism. We 
define the (electronic) Fock space as 

JV 



•^:=C®0/\L 2 (ft) 



N>1 1 

The grand canonical Hamiltonian is then given by Hq := Q) N>0 Hq with 
the convention that Hq = (1/2) X)i?^K'ez 3 nQ 1-^ — G C. The number 

operator reads A/" := © A r> N. It is then straightforward to check that 

E(Q) = inf a Fn (H n ) = inf tr ^ (JT n r) . 

rso^n), r =r, 

0<r<l, tr jr n (T)=l. 

The free energy at temperature 1//3 and chemical potential /i e M is defined 

by 

F(0,/3, M ):= inf (tr ^((if n - mAOH + ^tr ^(rio g r) 

o<r<i, tr ^ n (r)=i. 

= -Ilogtr^, [ e -«»n-/*AO] . (8) 

As explained in Introduction, our purpose is to prove that 

~|n Hoo e|fl| and F(fi, /?,/*) ~, nHoo /(/3, M )|£1| (9) 

in an appropriate sense. The first important property of E and F is the 
stability of matter. 
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Theorem 1.1 (Stability of Matter [16]). There exists a constant C 
such that the following holds: 

E(fl) > -C\Q\, F(Sl,0,n) > -C(l + /3- 5/2 +max(0,M) 5/2 ) |fi| 

for any bounded open set SlcR 3 and any {3 > 0, p e M. 

Sketch of the proof. The first step is to use an inequality for classi- 
cal systems due to Baxter [1], improved later by Lieb and Yau [28], and 
which allows to bound the full TV-body Coulomb potential by a one-body 
potential: 

v ( „,...,,„)>-£;3ews (10) 



i=l 




where S(x) = inf^ eZ 3 \x — R\ is the distance to the closest nucleus. Hence 
we have the lower bound 

_ 3/2 + V2 \ 
5(xi) J ' 

Next we split the kinetic energy in two parts and we use the uncertainty 
principle to show that on L 2 (Q) 

A 3/2 + V2 
4 6(x) ~ 

In proving this lower bound, one uses the Sobolev inequality in a small 
ball around each nucleus, exploiting the fact that the nuclei are fixed and 
separated by a distance at least one to each other. The proof of the stability 
of matter for systems with classical nuclei whose position is unknown is 
more difficult and it uses the improved version of (10) contained in the 
paper by Lieb and Yau [28], as explained in our work [16]. This shows 

<>E(-%- C ) hcncc Hn>--J2^-CM (11) 



i=i 



on L 2 (Q) and Tn respectively. The last step is to use the Lieb-Thirring 
inequality [27] which states that 

^(-A Xi )*,^ >C LT J^{xf'*dx (12) 

for all N > 1 and all Y-body fermionic wavefunction ^ S /\^ L 2 (Q). The 
density of charge is as usual defined by p^{x) = N J nN -i \^(x,y)\ 2 dy. 
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Using the fact that JqPv = N and Holder's inequality, (12) yields on the 
Fock space 

E(- A -) > c LT i^r 2/ w 5 / 3 . (is) 

Hence we obtain H n > (C LT /4)|0|- 2 / 3 7V 5/3 - CAf which, when optimized 
over N, gives the result for the ground state energy. 

For the free energy, we use (11), (13) and Peierls' inequality [30,34] to 

get 

F(p,»,Q) > -ilogtr^( e -^.(- A ^ 4 ) -C{l + $ 2 )\Sl\. 

The first term of the r.h.s. is the free energy of a free-electron gas which is 
bounded below by — C(l + /3 _5 / 2 )|f2| in the thermodynamic limit [16]. □ 

In order to state our main result, we need the 

Definition 1.1 (Regular sets in R 3 ). Let be a > and e > 0. We say 

that a bounded open set f2 C R 3 has an a-regular boundary in the sense of 
Fisher if, denoting by dSl = fi \ the boundary of ft, 

vte[o,i], g R 3 I d(x,dn) < \n\^ 3 t} < \n\at. (14) 

We say that a bounded open set SI C R 3 satisfies the £-cone property if for 
any x £ SI there is a unit vector a x G R 3 such that 

{y G K 3 | (x - y) ■ a x > (1 - e 2 )\x - y\, \x - y\ < e} C SI. 

We denote by !Z a .s the set of all SI C R 3 which have an a-regular boundary 
and such that both SI and R 3 \ SI satisfy the e-cone property. 

Note that any open convex set is in lZ a , e for some a > large enough 
and £ > small enough [15]. We may state our main 

Theorem 1.2 (Thermodynamic Limit for the Crystal [16]). There 
exist e <E R and a function f : (0, oo) x R — > R such that the follow- 
ing holds: for any sequence {Sl n } n >i C lZ a , e of domains with \Sl n \ — » oo, 
|r2n|~ 1/3 diam(£l„) < C, a > a > and < e < e 

E(Sl n ) _ F(n n ,/3,») - 
lim |o I =e, hm — =/(/?, m)- (15) 

Moreover f takes the form /(/3,/i) = <p(/3) — 
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Remark 1.1. We know from [22, Appendix A p. 385] and [10, Lemma 1] 
that if each set fi„ of the considered sequence is connected, then automat- 
ically |fi„|- 1 / 3 diam(f7„) < C. 

A very similar result was proved by C. Fefferman [9]. Our result is more 
general: we allow any sequence fl n tending to infinity and which is regular in 
the sense that {Q n }n>i Q 7Za,s- In Fefferman's paper [9], il n — £ n (fl + x n ) 
where £ n —* oo, il is a fixed convex open set and x n is any sequence in R 3 . 
These sets are always in lZ a . E for some a, e > 0. 

In our work [16] a result even more general than Theorem 1.2 is shown: 
we are able to prove the existence of the same thermodynamic limit if the 
crystal is locally perturbed (for instance finitely many nuclei are moved or 
their charge is changed) . A similar result can also be proved for the Hartree- 
Fock model. 



2. Other models 

Our approach [15,16] is general and it can be applied to a variety of models, 
not only the crystal case. We quickly mention two such examples. It is inter- 
esting to note that for these other models, we do not need the cone property 
and we can weaken the assumptions on the regularity of the boundary by 
replacing t on the r.h.s. of (14) by any t p , < p < 1. Details may be found 
in our article [16]. Roughly speaking, when the system is "rigid" like for 
the crystal (the nuclei are fixed), the proof is more complicated and more 
assumptions are needed on the sequence of domains to avoid undesirable 
boundary effects. 

2.1. Quantum particles in a periodic magnetic field. 

Define the magnetic kinetic energy T(A) = (— iW + A(x)) 2 where B = V x A 
is periodic (for instance constant) and A e L[ 2 oc (K 3 ). Next, consider the 
Hamiltonian 

N K 

H»> K := T(A) Xi + ]T T(A) Rk + V(x, R), 

i=l fe=l 

V(x, R) = Tr— —\ + 2 E \x-^r\ + 2 .5. Wk 



,_tfe — Xi\ L \Xi — Xj\ z — , \it k - Rk> I 



The ground state energy is this time defined as 
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We do not precise the symmetry S of the particles of charge z which can be 
bosons or fermions. A formula similar to (8) may be used for the free energy 
on the (electronic and nucleic) Fock space. We prove in our paper [16] a 
result similar to Theorem 1.2 for this model. Lieb and Lebowitz already 
proved it in the seminal paper [22] when A = 0. They used as an essential 
tool the rotation-invariance of the system to obtain screening. When A ^ 
the system is no more invariant by rotations and their method cannot be 
applied. 

2.2. Classical nuclei with optimized position. 

For all R C f2, #i? < oo, let us define 

Hg* + V{x,R) 

i=l 

and the associated ground state energy by 

#R<ic 

We could as well optimize the charges in [0, z] of the nuclei without changing 
the energy [6,16]. However, the free energy itself is not the same when the 
charges of the nuclei are optimized or not [16]. 

Surprisingly, to our knowledge the existence of the thermodynamic limit 
for this model was unknown. A result similar to Theorem 1.2 is proved in 
our paper [16] for E" . 

3. A general method 

In this section, we give the main ideas of our new approach which allows to 
prove Theorem 1.2 and its counterparts for the other models quoted before. 

3.1. Screening via the Graf-Schenker inequality 

As mentioned in the introduction, an important step is to quantify screen- 
ing. For quantum nuclei without a magnetic field (A = 0), Lieb and 
Lebowitz used [22] the following method (see Figure 1). First they took a 
big ball B which they packed with several small balls of different size. In 
each of these balls, they took the (neutral) ground state of the correspond- 
ing ball. As the system is invariant under rotations, they can freely rotate 
each ground state. Averaging over rotations of all the small balls, they re- 
duced the computation of the interaction between them to that of classical 
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Fig. 1. A comparison between the original method of Lieb and Lebowitz [22] (left) and 
our method based on the Graf-Schenker inequality [13,15,16] (right). 



pointwise particles located at the center of the balls, by Newton's theorem. 
As each subsystem is neutral, this interaction vanishes. This proves a for- 
tiori that there exists an adequate rotation of each system in each little 
ball such that the total interaction between them cancels. Choosing this 
configuration, they could build a test function whose energy is just the sum 
of the small energies, proving an estimate of the form E(B) < ^ fe E(Bk). 
This inequality can be used to prove the limit for balls. Clearly this trick 
can only be used for rotation-invariant systems. 

Note in the Lieb-Lebowitz proof, a domain (the big ball) is split in 
several fixed subdomains and an average is done over rotations of the states 
in each small domain. This yields an upper bound to the energy. The Graf- 
Schcnkcr inequality is kind of dual to the above method (see Figure 1). 
This time a domain £1 is split in several subdomains by using a tiling of the 
space R 3 . But the system is frozen in the state of the big domain and the 
average is done over the position of the tiling. This yields a lower bound to 
the energy of the form E(Q) > J2k E(A^ r ' u ^ fl Q) + errors, where A[ r,u ^ are 
the tetrahedrons which make up the (translated and rotated) tiling. 

The Graf-Schenker inequality was inspired by previous works of Conlon, 
Lieb and Yau [4,5]. It is an estimate on the Coulomb energy of classical 
particles. The proof of Fefferman in the crystal case [9] was also based on a 
lower bound on the free energy in a big set and an average over translations 
of a covering of this set (the method was reexplained later in details by 
Hugues [17]). Fefferman [9] uses a covering with balls and cubes of different 
size. The lower bound depends on the number of balls contained in the big 
domain and of the form of the kinetic energy which is used to control error 
terms. 
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Let G = R 3 x SO3 (R) be the group of translations and rotations acting 
on R 3 , and denote by dX(g) its Haar measure. 

Lemma 3.1 (Graf-Schenker inequality [13]). Let A be a simplex in 
M 3 . There exists a constant C such that for any N G N 7 Z\,...,Zn € R, 
Xi G R 3 and any £ > 0, 

\ - ZjZj f d\(g) \ - ZiZjtg£A(xi)t g eA(xj) C h 



Xi X j 



J G \£A\ ^ \xi-xA i^ Zl ' 



l<i<j<N '~* " Jl "° 1 l<i<j<N 1 ' " i=l 

(16) 

In the previous theorem it is not assumed that A yields a tiling of R 3 . 
Up to an error which scales like £, (16) says that the total Coulomb energy 
can be bounded from below by the Coulomb energy (per unit volume) 
of the particles which are in the (dilated) simplex g£A, averaged over all 
translations and rotations g of this simplex. 

Because of the above inequality simpliccs play a specific role in the study 
of Coulomb systems. Hence proving the existence of the thermodynamic 
limit for simplices first is natural (as it was natural to consider balls in the 
Lieb-Lebowitz case due to the invariance by rotation). In the next section 
we give an abstract setting for proving the existence of the limit when an 
inequality of the form (16) holds true. 



3.2. An abstract result 

In this section we consider an abstract energy E : CI S A4 i— » E(Cl) e R 
defined on the set M of all bounded open subsets of R 3 and we give sufficient 
conditions for the existence of the thermodynamic limit. In the application, 
E will be either the ground state energy or the free energy of the system 
under consideration. 

We fix a reference set A e TZ a ,e which is only assumed to be a bounded 
open convex set in R 3 (it need not be a simplex for this section), such that 
e A. Here a,e > are fixed. We assume that the energy E satisfies the 
following five assumptions: 

(Al) (Normalization). E(<fr) = 0. 

(A2) (Stability). VO e M, E(Sl) > -n\fl\. 

(A3) (Translation Invariance). VO € TZ a .s, Vz G Z 3 , E(fl + z) = E(Ct). 
(A4) (Continuity). Vf2 e ft a , e , ft' e 1l a >^ with ft' C ft and d(dfl, dfl') > S, 

E(si) < E(n') + K \n \ n'\ + \n\a(\n\). 
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(A5) (Subaverage Property). For all ft € A4, we have 

E(Q) > £ £?(n n .g • (*A)) d\(g) - |fi| r a(*) (17) 

where |fi| r := inf{|0|, f2 C f2, Sle 7\L a , £ } is the regularized volume of SI. 

In the assumptions above a is a fixed function which tends to at 
infinity and S,a',e' are fixed positive constants. In our work [15], an even 
more general setting is provided. First (A3) can be replaced by a much 
weaker assumption but we do not detail this here. Also a generic class of 
regular sets 1Z is considered instead of TZ a ,e- This is because for instance 
the cone property is only needed for the crystal case and it is not at all 
necessary in other models, hence the concept of regularity depends on the 
application. 

Notice (A4) essentially says that a small decrease of SI will not deacrease 
too much the energy. A similar property was used and proved in the crystal 
case by Fefferman [9, Lemma 2]. Taking SI' = and using (Al), property 
(A4) in particular implies that for any regular set SI <G lZ a , 6 , E(SY) < C\Sl\. 
However this upper bound need not be true for all SI G M . We give a sketch 
of the proof of the following result in Section 3.5. 

Theorem 3.1 (Abstract Thermodynamic Limit for A [15]). 

Assume E : M — > K satisfies the above properties (A1)-(A5) for some 
open convex set A € TZ a ,e with € A. There exists e 6 R such that ee(g) = 
\£A\~ 1 E(g£A^j converges uniformly towards e for g e G = K 3 x SO(3) and 
as £ — > oo. Additionally, the limit e does not depend on the set A*. 

3.3. Idea of the proof of (A1)-(A5) for the crystal 

Before switching to the abstract case of a general sequence {0„}, we give 
an idea of the proof of (A1)-(A5) in the crystal case. We apply the theory 
of the previous section to both the ground state energy and the free energy 
of the crystal which were defined in Section 1. First (Al) and (A3) are 
obvious. Property (A2) is the stability of matter as stated in Theorem 1.1. 
On the other hand (A5) is essentially the Graf-Schenker inequality (16), 
up to some localization issues of the kinetic energy which essentially have 
already been delt with by Graf and Schenker [13]. 

For the crystal the most difficult property is ( A4) . The difficulty arises 
from the fact that this is a very rigid system. For the two other examples 



*This means if all the assumptions are true for another set A' then one must have e' = e. 
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Fig. 2. Idea of the proof of (A4) for the crystal. 



mentioned in Section 2, (A4) is obvious, the energy being nonincreasing: 
E(Cl) < E(W). This is because we can simply choose a ground state of il' 
as a test for f2 and take the vacuum in £1 \ Q'. In the crystal case we always 
have nuclei in Q \ tt' and if we do not put any electron to screen them, they 
will create an enormous electrostatic energy. 

The idea of the proof of (A4) for the crystal is displayed in Figure 2. We 
build a test state in f2 by considering the ground state in f2', and placing 
one radial electron in a ball of fixed size on top of each nucleus ouside f2'. 
By Newton's theorem, the electrostatic potential out of the support of the 
electron will vanish, hence the energy will simply be E(Q') plus the sum of 
the kinetic energies of the electrons, which is bounded above by a constant 
times |f2 \ for regular domains. The only problem is that we cannot put 
an electron on top of the nuclei which are too close to the boundary of fl 
or of Q' . For these nuclei, using the cone property we can place the ball 
aside and create a dipole. The difficult task is then to compute a bound on 
the total interaction between the dipoles and the ground state in £1' . We 
prove [16] that it is o(|f2|), using a specific version of stability of matter. 



3.4. General domains and strong subadditivity of entropy 

In the previous two subsections, we have presented our abstract theory giv- 
ing the thermodynamic limit of special sequences built upon the reference 
set A, and we have explained how to apply it to the crystal case. For all 
regular domain sequences we can only get from (A5) a bound of the form 

. E(n n ) _ 

hmmf > e. 

ri^oo | Sin 

In order to get the upper bound, we use a big simplex L„A of the same size 
as Q n and a tiling made with simplices of size £ n <C L n , as shown in Figure 
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3. We use the ground state of the big simplex L n A to build a test state in 
f2„, hence giving the appropriate upper bound. To this end, we need some 
localization features, hence more assumptions in the general theory. 




Fig. 3. Proof for general sequences {f2 n }. 



It is sufficient [15] to assume that 

(i) A can be used to build a tiling of K 3 ; 

(ii) the free energy is essentially "two-body" t such that we may write the 
total energy E(L n A) as the sum of the energies of the small sets of the 
tiling, plus the interaction between them and the relative entropy; 

(Hi) the entropy is strongly subadditive. 

This is summarized in the following assumption. We assume that T is a 
subgroup of G yielding a tiling of R 3 by means of A, i.e. U^gr/^A = R 3 
and //A n vA = for /i ^ v. 

(A6) (Two-body decomposition). For all L and I we can find g € G and 
maps E g : T -> R, I g : T x V -> R, s g : {V : V C T} -> R such that 

• £ ff 0) = / ff 0, i/) = if %iA n (LA) = 0; 

. E{LA)>Y,EM + \ E y^)-s 9 (r)-|£A|aW; 

• For all C T and A v = LA n U^ep A 



t We could as well assume that the energy is fc-body with A; < oo but this would complicate 
the assumptions further more. 
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• (Strong subadditivity). for any disjoint subsets Pi, P 2 , V3 C T 



In the applications* the previous quantities are interpreted as follows: 
E g (V) is the free energy in the union Ap = (LA) n U^-pig^A, I g (n, v) is 
the interaction energy between the simplices lg[iA and IgvA, and s g (V) is 
the difference between the entropy of Ap and the sum of the entropies of 
igfiA with fieV. 

Conjectured by Lanford and Robinson [18] the strong subadditivity 
(SSA) of the entropy in the quantum mechanical case was proved by Lieb 
and Ruskai [23,24]. The fact that SSA is very important in the thermody- 
namic limit was remarked by Robinson and Ruelle [31] and others [34]. In 
the article [15] we prove the following 

Theorem 3.2 (Abstract Limit for general domains [15]). Assume 
E : M. — > M satisfies the properties (A1)-(A6) for some open convex 
polyhedron A G lZ a , e with G A, yielding a tiling 0/R 3 . Then we have for 
all sequences {fl n } C TZ a ,E with \fl n \ — ► 00 and |Q n | _1 / 3 diam(fi„) < C, 



where e is the limit obtained in Theorem 3.1. 

The proof of Theorem 3.2 is based on a careful estimate of the energy 
and the interaction energies of boundary terms, ie. of the sets lg\xA which 
intersect the boundary of the big set LA. The application to the crystal 
is not much more difficult than for Theorem 3.1. Indeed in the paper of 
Graf and Schenker [13], (16) was expressed using a tiling of M 3 and the 
last subaverage property of (A6) essentially follows from their ideas [13]. 
Strong subadditivity of the entropy is usually expressed via partial traces. 
A generalization in the setting of localization in Fock space is detailed in 
our article [16]. 



s g (Pi UP 2 U Vz) + s g (V 2 ) < s g (Vx U V 2 ) + s g (V 2 U V 3 ); 




G/r 




= e 



tDuc to some localization issues of the kinetic energy, it is often needed that the sets of 
the tiling slightly overlap. See [15] for a generalization in this direction. 
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3.5. Proof of Theorem 3.1 

Denote as in the Theorem e e (g) = E(g£A)\£A\~ 1 . Notice that (A2), (A4) 
with O' = 0, and (Al) imply that eg is uniformly bounded on G. Also we 
have by (A3) e e (u + z,R) = e e (u, R) for all (u, R) G R 3 x SO s (R), z e Z 3 , 
i.e. is periodic with respect to translations. Hence it suffices to prove the 
theorem for g = (u, R) E [0, l] 3 x S0 3 (R). 

Next we take g e G, L » £ and apply (A5) with f2 = gLA. We get 

, , l-a(£) f E(gLAC\ g£A) , 

eLi§) * -JlAT J g \£A\ ds ~ ^ 

Let us introduce the set Z of points z e Z 3 such that i?€A + u + z C g£A 
for all u £ [0, l] 3 and all R € SC^R). We also define <9Z as the set of points 
z e Z 3 such that (R£A + u + z)DgLA ^ for some (u, R) E [0, l] 3 x S0 3 (R) 
but z ^ Z. We obtain using (Al) and (A3) 

E(gLAf\g£A) ^ f f E(gLA n (MA + u + z)) 



rE {9 LAng£A) dg= r d r dR 

Jg l« A l zl^-Vi] 3 Jso 3 (r) 



|*A| 



V /' dti /' 



.^(ffLAnfWA + n + z)) 
art 



+(#Z) [ du I dR e e (u,R). 

J\0,1] 3 JSOitK) 



I [0,1] 3 JS0 3 (R) 

Using the stability property (A2), we infer 

E(gLA n (MA + u + z)) \gLA n (MA + u + z)\ 

i^Aj - ^ j/Aj - K - 

Hence 

/ ^ A n , #) ^>(P)f e,( 3) ^ + K (#aZ). 

Jg J[04] 3 xso 3 (r) 

As A has an a-regular boundary, it can be seen that (#dZ) < CL 2 £ and 
#Z = \LA\+0(L 2 £). Using again that ep is bounded, we eventually obtain 
the estimate 

e L (g)> f e e (g)dg-C(a(£)+£/L) 

J[0,l] 3 xSO 3 (R) 

for some constant G. It is then an easy exercise to prove that 
lim inf eg = lim 

<->oo G l ^°°J[0,l] 3 xSO 3 

and finally that eg — > e in i 1 ([0, l] 3 x SO s 
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The last step consists in proving the uniform convergence, using (A4). 
Fix some small rj > 0. As £ A and A is convex, we have (1 — 77) A C A. 
More precisely, there exists an r > and a neighborhood W of the identity 
in S0 3 (R) such that R(l-r])A + u C A for all (u,R) £ A := B(0,r)xW C 
G. We have that g£(l - 77) A C £A for all g £ A t := B(0,r£) x W, hence 
in particular for all g £ A. Now we fix some g £ G and apply (A4) with 
ft = glA and ft' = <?^(1 - 77) A, we get 

E(g£A) < E{gg£{\ - rj)A) + C\£A\rj + o(\£A\). 

Integrating over g £ A and dividing by \£A\ we infer 

ee(g) < J ^ e{i- v )t{g) dg + Cr 1 + 0(1)^00. 

First we pass to the limit as £ — > 00 using that — > e in i 1 (G) and |A| 7^ 0. 
Then we take 77 — > and get limsup^^^ supg gG eg(g) < e. This ends the 
proof of Theorem 3.1. □ 
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